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Abstract 

The cubic Henon-Heiles system contains parameters, for most values of which, the 
system is not integrable. In such parameter regimes, the general solution is expressible 
in formal expansions about arbitrary movable branch points, the so-called psi-series 
expansions. In this paper, the convergence of known, as well as new, psi-series solutions 
on real time intervals is proved, thereby establishing that the formal solutions are 
actual solutions. 

1. Introduction 

The generalized Henon-Heiles Hamiltonian is 

H = - (x 2 + y 2 + Ax 2 + By 2 ) + Dx 2 y m - 2 - -y m , (1.1) 
2 m 

where x = t?, y = M and A, B, C and D are nonzero parameters. It corresponds to 
the original form given by Henon and Heiles [1] in their studies of astronomical problems 
when A = B = C = D = 1 and m = 3. Here, we consider the cubic form (m = 3) of (1.1), 
with the corresponding equations of motion 

x + Ax + 2Dxy = 0, 

ij + By + Dx 2 - Cy 2 = 0. 

Several studies have been made regarding the integrability properties and singularity struc- 
tures of various versions of (1.2); for example, Chang, Tabor and Weiss [2] have analy- 
sed (1.2) with A = B = 1, whereas Chang, Greene, Tabor and Weiss [3] have considered 
the same system for general m. In keeping with their notation, we rescale the variables, 

x y \ D n q\ 
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to express (1.2) as 

x + Ax + 2Xxy = 0, 

y + 5y + Xx 2 - y 2 = 0. 

The leading-order behaviours of solutions about an arbitrary singularity, as well as their 
resonance structures (in the sense of the ARS algorithm, Ablowitz, Ramani and Segur [4, 
5]), have been derived (see [2, 3]; Chang, Tabor, Weiss and Corliss [6]; Bountis, Segur and 
Vivaldi [7]; Weiss [8]; Tabor [9, p. 337]). 

We summarize the results briefly, in order to fix notation and to single out the present 
objectives. Setting r = t — to, where to is an arbitrary constant (the location of the 
singularity), and 

x ~ ar a , y ~ br?, t -► 0, (1.5) 
two types of leading-order singular behaviours are found: 

(i) a = /3 = -2, 

(1.6 

(ii) (3 = -2, Re(a) > (3. y ' 
Carrying the analysis further by writing 

x ~ ar a +pr a+kr , y ~ br a + qr p+kr , (1.7) 

and determining those values of k r for which p or q is undetermined, one finds the corre- 
sponding resonances. It turns out that 

(i) a = (3 = — 2, k r = — 1, 6, r, f, 



5 1/ 1 
where r,f = - ± -\1 - 24(- + 1); 

2 2V (1.8) 
11 

(ii) Re(a)>-2, (3 = -2, a,a = - ± - 48A, A; r = -1, 6, 0, r, f, 

where r, f = =F\/1 — 48A. 

In case (ii), a and a are two possible leading orders, with corresponding resonances r 
and f. This case occurs only for A > — |, since i?e(a) > — 2. 

The nature of the resonances for each value of A is summarized as follows: 

Case(i) 

24 

A < — — : r and r are complex 

24 5 

A = : r = r= - 

23 2 

24 1 

— < A < -- : r>0, r>0, r^r 



23 2 
1 

~ ~2 



r = 5, r = 



--<A<0 :r>0, f < 
2 

A > : r and f are complex 
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Case (ii) 

f > 
f > 


f are pure imaginary 

Those cases where a negative resonance is present (in addition to —1) likely correspond 
to singular, rather than general, solutions (see, e.g., [9, p. 339]). The significance of 
repeated resonances (case (i), A = — §3, and case (ii), A = ^) is, at present, unknown. 
The case A = — \ is somewhat anomalous [9, p. 340] and we shall comment upon it later. 
In those cases where there exist four distinct resonances, the usual —1 and three others 
with nonnegative real parts, general solutions may be written down in series expansions 
about the arbitrary singularity to- These may be Laurent series, if the system passes the 
Painleve test [5] or, otherwise, psi-series (Hille [10, p. 249]). 

It has been found [2] that the system is integrable when A = — 1 and A = B, when 
A = — jjr, and when A = — jq and B = 16 A. In the latter case, although the system is 
integrable, it passes only the "weak" Painleve test in that it admits a rational movable 
branch point. Outside of these three cases, the system possesses solutions with movable 
branch points (logarithmic, irrational or complex), and a psi-series expansion is required to 
represent the general solution about an arbitrary branch point. Such expansions constitute 
formal general solutions containing four arbitrary constants, the remaining coefficients 
being well defined by a self-consistent recursion relation. 

In this paper, we prove the absolute convergence of such series on real intervals of the 
form < r < R, R > 0, thus establishing that the formal solutions are actual solutions. 
Similar results have been obtained by Melkonian and Zypchen [11] for the Lorenz system 
(see, e.g., [9, p. 344] and Sparrow [12]), where only logarithmic psi-series occur. A study 
of the convergence problem for a class of second- and third-order ordinary differential 
equations has been made by Hemmi and Melkonian [13]. Some of the key ideas within 
the proofs may be extracted from a paper by Hille [14] regarding second-order quadratic 
systems. Studies of general systems of partial differential equations which admit WTC 
(Weiss, Tabor and Carnevale [15]) expansions containing logarithms have been made by 
Kichenassamy and Srinivan [16] and Kichenassamy and Littman [17], who have proved 
convergence by entirely different methods. Logarithmic psi-series occur for the Cubic 
Henon Heiles system (1.4) when A = — 1 and A / B, or when A = — jq and B / 16 A, and 
these may be dealt with by the methods discussed in [11, 13, 16, 17]. 

Here, we restrict our attention to the non-logarithmic cases. Sections 2 and 3 deal 
with the case(i)-leading orders, involving a series with complex exponents (Section 2) 
and one with irrational exponents (Section 3). The necessary Lemmas are relegated to 
Appendix A. The series which we employ in Sections 2 and 3 are not the same as the 
ones given, e.g., in [3], and we confirm the validity (self-consistency) of our version in 
Appendix B. Section 4 concerns the case(ii)-leading orders where, once again, we employ 



1 

— < A < 
2 


: r 


<o, 


< A < : 
48 


, r 


<o, 




. r 




»>5 


. r 


and 
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a new series, the validity of which is confirmed in Appendix C. Concluding remarks are 
made in Section 5, including a discussion of the case A = — \. 



2. Case(i)-leading orders: ex. = (3 = —2, 



k r = -1, 6, | ± \Jl - 24 (\ + 1) complex 



Here, the resonances r, r = \ ± 1 — 24(^ + 1) are complex, with A < — 1| or A > 0. 
The general solution may be taken in the form 



oo oo oo 

k— 2+rl+fm 
m 1 ; 



= E E E akw ' 

(2.1) 



x 

k=0 1=0 m=0 

oo oo oo 

y = YYY b ^ k ' 2+rl+rm - 



k=0 1=0 m=0 

The self-consistency of (2.1), as well as the presence of four arbitrary constants, is shown 
in Appendix B. It is also shown therein that (2.1) is equivalent to 



oo oo 

k-2+fl 



= YY [ akirk 2+rl + ^ kiT 

k=0 1=0 

oo oo 



(2.2) 

k-2+rl 



X 

k=0 1=0 

oo oo 

y = 

k=0 1=0 

and that (2.2) is equivalent to the solution given in [3]. Either (2.1) or (2.2) may be 
employed to prove convergence. However, the advantage of starting with (2.1) is that it 
is much easier to deal with a single, triply-indexed series such as (2.1) than a sum of two, 
doubly-indexed ones such as (2.2), vis-a-vis the confirmation of its validity. 

The convergence proof consists of resumming (2.2) into more tractable forms, and 
showing that the latter are majorized by series which converge on an interval, by the ratio 

test. 



. Thus, let e = \^2A{\ + 1) - 1 > 0, note that r - 1, r - 1 = | ± ie, and write 

oo oo 

EE [akir k - 2+l HI^) 1 + a kl T k - 2+l+ (l- l£ y 



oo oo 

X = 

::=0 1=0 

oo oo 



(2.3) 



E E t( ajW + afe ^ cos ( £lz ) + Kaki - an) sin(dz)] e^ lz T h 



k-2+l 

fc=0 1=0 

where z = hi(r), r > 0. It follows that 



x = ^f y (z)r^, (2.4) 

7=0 

where 

f~f( z ) = E K a ki + a k i) cos(elz) + i(a ki - a M ) sm(elz)} e^ lz (2.5) 

fc+/=7 
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is a polynomial in e^ z of degree 7, with coefficients that are bounded functions of z. 
Similarly, 

00 

y = Y J h 1 {z)r^\ (2.6) 

7=0 

where 

M z ) = i( bkl + bki) cos(elz) + i(b kl - b M ) sm(elz)] ei lz (2.7) 

fe+/=7 

is a polynomial of the same type as / 7 (z). 

Let u(t) = x(t) = % and v(t) = y(t) = §, then 

00 00 

u = ± = x>(*K 7 - 3 = E [(7 - 2)/ 7 + /;] r<-\ 
7=0 7=0 

00 00 

v = y = Y J k,(z)r^ = [(7 - 2)^ + ^] ^~ 3 , 
7=0 7=0 



(2.8) 



where / 7 = ^,/i 7 = and g 7 and /c 7 are also polynomials of the same type as /, 



7 



and hj. 



Express (1.4) as a system of four first-order equations in x,u,y,v to obtain 

(2.9) 



x = u, 

u = —Ax — 2Xxy, 
y = v, 
v = —By — Xx 2 + y 2 . 



Substitution of (2.4), (2.6) and (2.8) into (2.9) gives, for 7 > 0, 
/; + ( 7 -2)/ 7 -<7 7 = 0, 

7 

<? 7 + (7 - 3) 57 + A/ 7 _ 2 + 2A fi-i*K = °> 

+ (7 - 2)/l 7 - fey = 0, 

7 7 

fc 7 + (7 - 3)/c 7 + B/i 7 _ 2 + A ^2 f^-nfn - ^2 h i-n h v = °- 
For 7 = 0, we find 



(2.10) 



3/1 36 

/ = ±-y- + 2, 50 = -2/o, ^0 = -^, k = j. (2.11) 

For 7 > 0, (2.10) is expressible as 

4 + A 7 /; = F 7 , (2.12) 
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where 



\ fc 7 / 



/ 7-2 -1 0\ 
-6 7-3 2A/ 
7-2 -1 



V 







2A/ 



A 



7-3 



V 



7-1 

-A/ 7 _ 2 -2A^/ 7 _ M ^ 



7— 1 7— 1 

-Bhy-2 — A ^ fj-nfn + ^ h-j-^h^ 

jU=l /1=1 



(2.13) 



/ ^7 \ 

v * 7 y 



The eigenvalues of A 7 are 7 + 1,7 — 6 and 7 — (| ± ie), precisely 7 — &v where /c r is a 
resonance. The matrices A 7 are simultaneously diagonalizable by a matrix P independent 
of 7. -P is the matrix of eigenvectors, 



/ -A/o -A/o 
3A/ -4A/ 



A/o 



V 



A/o ( g + * e 



3 2 + 



1\ A 



A/o 



A/o 
1 



12 3 2 + - - + is) 3 2 + - 



3 2 + 



1\ I 



is 

2 



(2.14) 



and P l A^P = D 7 , where 

/ 7+ 1 

7-6 



V 








7 









+ IE 



7-1 



(2.15) 



7 



For 7 > 6 (in which case all eigenvalues of A 7 have positive real parts), the solution of 
the matrix differential equation (2.12) may be expressed as 

f 7 (z) = f* e A ^ x - z ^F y {x)dx = f* Pe D ^ x - z ^p- 1 F 1 {x)dx, (2.16) 
J — 00 J — 00 

where the exponential within the integrals denotes the exponential of a matrix. 
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Given an m x n matrix A = (ay), define its norm as 

\A\\ = max I } ' \aij\ I . (2.17) 




It follows from (2.16) that 

||/;^)|| < IIPIHIP- 1 !! f e^- 6 )^)||F 7 (x)||dx, (2.18) 

J — CO 

since \\ e D-,(x-z) y = e { 7 -6)(x-z) for x < z . p depends upon A, so that ||P||||P _1 || < M (A), 
a constant (independent of 7). Thus, 

< M(A) ( Z e^- 6 ^ x -^\\F^{x)\\dx. (2.19) 



Theorem 2.1. T/iere exisfe K > suc/i i/tai /or a// z < and 7 > 1, 

< -^2=. (2.20) 

3 

Proof. Since / 7 (z) is a polynomial in X = ea 2 of degree 7, with coefficients that are 
bounded functions of z (linear combinations of cos(elz) and sin(e/z)), 

|/ 7 (*)|< J>We!™ = P 7 (X), (2.21) 

m=0 

with Om > for < m < 7. By Lemma Al of Appendix A (with q = 1, ra 7 = 7, M = 1 
and p = 1), given A > 1, there exists > such that 

|/ 7 (*)| < P 7 (X) < {Kl J^ V for 1< 7 <A-1. (2.22) 

Similarly, there exist positive constants K2, A3 and K4 corresponding to the polynomials 
g-y, /i 7 and /c 7 , repectively. With K = max{ATj}, we obtain 

l<i<4 

||/ 7 (z)|| < L_^_i_ < for 1< 7 <A-1, z<0, (2.23) 

3 

since < e^ z < 1. The inequality (2.23) constitutes our inductive hypothesis, and we 
shall show that (2.20) holds for 7 = N, thus establishing that it holds for all A > 1. 
By (2.13) and the inductive hypothesis (2.23), we find that 

N-l 

\G N \ <\A\\f N „ 2 \+2\\\Y,\fN^\\K\ 

(U=l 
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|A V | <V B \t K) * +(\M + 1)*\(2K) N . (2.25) 



by Lemma A2 of Appendix A. Similarly, 

B\mr 

Thus, 

\\F N \\ < E(2K) N , (2.26) 

where 

E < max{\A\ + 2|A|vr, \B\ + (|A| + 1)tt}. (2.27) 
Thus, from (2.19) and (2.26), we obtain 



\\fN(z)\\ < M(X)E l* e^ N - e)(x - z) (2K) N dx 

J — oo 



M(X)EVNTT\ (2Kf for 



A-6 J v / ATl 



(2.28) 



For all sufficiently large N, MiQE^E EL < x, so beginning with such an N > 6, we obtain 

IIM*)II < ±=L=, (2.29) 

which completes the proof. 
By (2.4), 

oo 

x = ^2M z y- 2 , (2-30) 

7=0 

and by (2.20), 

\m^ 2 \ < |=^- 2 , (2.3i) 

and the series 

(, 32) 

converges absolutely by the ratio test for r < ^ . Thus, the series (2.30) for x converges 
absolutely for < r < R, where R is at least ^ . Similarly, the series for x, y and y 
converge absolutely for < r < R. 
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3. Case(i)-leading orders: ol = (3 = —2. 



k r = -1, 6, | ± |y 1 - 24 (i + 1) irrational 

In the present case, the resonances r,f = | ± i^/l — 24 + l) are irrational, positive 

and distinct, with — 1| < A < — |. The general solution is given by (2.2) (or (2.1)) as in 
the complex-resonance case of Section 2, but the details of the resummation are slightly 
different since, in the present case, r < 1 and/or f < 1 is possible. 

Let n be the least positive integer such that fi\ = r — ^ > and ^2 = r — ^ > 0. (Any 
positive integer satisfying the latter two conditions will also do.) Then the solution x in 
(2.2) may be expressed as 



00 00 

x = E E h^~ 2+ " W + awr*-^^ 1 ' 

fc=0 1=0 

Let w; = r 1 /™ and z = ln(r), r > 0, then 

00 00 

*=££[< 



a M w nk - 2n+l e^ lz + auw^-^e^ 



k=0 1=0 

00 



= £ £ 

7=0 nk+l=j 



a kl e^ lz + a kl e^ lz ) ^- 2 " = ]T/ 7 ( 

7=0 



7-2n 



where 



f J (z)= £ (a kl e^ + a kl e^ lz ) 



nfc+i=7 

is a polynomial of degree 7 in the two variables e^ iz and e^ 2Z . Similarly, 



(3.1) 



(3.2) 



(3.3) 



7=0 



where 



is a polynomial of the same type as fj(z). 
Let 



u = x 



00 


00 


7=0 


7=0 


00 

EM^ 7 ~ 3n 

7=0 


00 

-£[(* 

7=0 



w 



7— 3n 



7— 3n 



(3.4) 



(3.5) 



(3.6) 
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where, as before, x = fj = ^j, etc., and g 1 and /c 7 are polynomials of the same type 
as / 7 and h 1 . Substitution of (3.2), (3.4) and (3.6) into the system (2.9) gives, for 7 = 0, 



^ = ± aVX + 2, So = -2/o, h o = -j, 
(as (2.11)) and, for 7 > 0, 

fj + ^-7/7 -^7 
(as (2.12)), where F 7 is again given by (2.13), but now, 

\ 



k = - y 



(1-2 -1 

n 



-6 1-3 2A/ 

n 







V 





2A/ 








n 



- - -3 

A n 



(3.7) 



(3.8) 



(3.9) 



and coincides with the A 1 of (2.13) if n = 1. The eigenvalues of A 7 in the present case 
are - — k r , where k r = — 1, 6, r, r are the resonances. 

The matrix P which diagonalizes all of the A 7 is the same as (2.14) which, for the 
present situation, is expressed as 



/ -A/o -A/o 
3A/ -4A/ 



P = 



3 

-9 



3 
12 



with P _1 A 7 P = D 7 and 



D 1 = 



( 1 + 1 

n 









7 
n 



\ 








A/o 
A/o(r-2) 

3(2 + 1 



A/o 
A/o(f - 2) 

3(2 + 1 



\ 



3(2 + -)(r-2) 3(2 + T )(f-2) 








7 

r 

n 



\ 










7 



— r 



n 



(3.10) 



(3.11) 



As before, ||P|| \\P 1 || is bounded by a constant M(A), the solution of the system (3.8) is 
expressible as 



/z 
p e Dt(*-*)p-iFi(x)da 
-00 



(3.12) 
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for ^ > 6 (noting that within the range of A presently under consideration, < r, r < 6, 
so that all eigenvalues of A 7 are positive for ^ > 6), and 

< M(A) f e(^ 6 )(— >||F 7 (s)||cir. (3.13) 

J — oo 

Theorem 3.1. T/iere exisfe K > suc/i i/iai /or a// z < and 7 > 1, 

H£(*)ll * J=Sf < 3 - 14 ) 

Proof. Since / 7 (z) is a polynomial of degree 7 in X = e MlZ and 1" = e M2Z , Lemma A3 of 
Appendix A (with M = 1 and p = 1) implies that given N > 1, there exists Ki > such 
that 

\f,(z)\< [ 1+ fOT !<7<iV-l. (3.15) 

Similarly, there exist positive constants K2, K3 and K4 corresponding to the polynomials 
g 7 , /i 7 and /c 7 , respectively. With i^T = max^{i^j}, we obtain 

\\m\\< { ^ 7 + 1 fOT !<7<iV-l. (3.16) 

Since < e^' 2 < 1 for i = 1,2 and for all z < 0, we have 

< for 1< 7 <AT-1. (3.17) 

The inequality (3.17) constitutes our inductive hypothesis, and we shall show that (3.14) 
holds for 7 = N, thereby establishing that it holds for all N > 1. 
As in Section 2, we find that 

II^tvII < E(3K) N , (3.18) 

where E is bounded as in (2.27). Thus, from (3.13) and (3.18), we obtain 

\\f N {z)\\ < M(X)E f Z e(^~ 6 ) ix - z) {3K) N dx 
J —00 



M(X)EVN + T \ (3K) n for N >6 
f - 6 J v^VTT n > 



(3.19) 



For all sufficiently large iV (and any n > 1), Mwj^m < i, so be 

ginning with such an 

n 

N > 6n, we obtain 

||/*(*)|| < (3.20) 

which completes the proof. 
By (3.2), 

00 

x = ^M^W~ 2n , (3-21) 
7=0 
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and by (3.14), 



\MzW- 2n \ < -pH=^- 2 ™, (3.22) 



and the series 
(MO 



y_§£L w 7-2n (3 . 23) 



7=0 



converges absolutely by the ratio test for w < j^, i.e., for r < ^ 3 ^„ ■ Thus, the series 

(3.21) for x converges absolutely for < r < R, where R is at least j^yi ■ Similarly, the 
series for x, y and y converge absolutely for < r < R. 

4. Case(ii)-leading orders: Re(a.) > —2, (3 = —2 

As mentioned in Section 1, if Re(a) > —2 and A > — |, two types of leading-order and 
resonance structures are possible, namely, 

(a) Leading order a = - + -y/l — 48A , resonances k r = — 1,6, 0, r; 

(b) Leading order a = - — — 48A , resonances /c r = — 1,6, 0,r, 

where r, f = =F\/1 — 48A. For A < r < 0, so that (a) does not correspond to a general 
solution, and will therefore not be considered. (See Conte, Fordy and Pickering [18] for 
a discussion of negative resonances and the appropriate series expansions.) In case (b), 
f>0for— i<A<0 and < A < ^, and f is pure imaginary for A > ^. The general 
solution takes the form 



oo oo oo 



fc=0 1=0 m=0 

4.2 

oo oo oo 
fc=0 1=0 m=0 

discussion of the validity of which, as well as other relevant issues, is relegated to Ap- 
pendix C. If f is pure imaginary (A > ^), then the methods employed earlier do not 
apply to (4.2) (see Section 5), and we have not determined whether the series converge or 
not. For f > (A < ^), we proceed as in the preceding sections. 

Let n be the least positive integer such that (3 = f — ^ > 0, let w = T X l n and z = 
ln(r), r > 0, in order to express x in (4.2) as 

OO OO oo 

k+2m+H+a+am+f3l 



X = > > > aklmT 
k=0 1=0 m=0 



E E E aklr ' 

(4.3) 

^ oo 
EE E a k lmW nk+2nm+l+ " n T° m+ P l = Y^f 1 {z)w^ n , 



k=0 1=0 m=0 7=0 
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where 



amz+f3lz 



fj( z ) = a klv,( 
nk+2nm+l=f 



is a polynomial of degree 7 in the two variables e az and e@ z . Similarly, 



(4.4) 



7=0 



where 



amz+fjlz 



nk+2nm+l=f 



is a polynomial of the same type as f-y(z). 
Let 



u = x 



00 00 

7+an- 



7=0 

00 



7=0 



v = y = Y, M*)*" 7-3 " = 2) ^ + K 



w 



w 



( 7-3n 



7+an-rt 



7=0 7=0 

and substitute (4.3), (4.5) and (4.7) into the system (2.9) to obtain, for 7 = 0, 

f arbitrary, g = af , h = 6, k = -12, 
and for 7 > 0, the system 

where, as before, 



h = 



and here, 



V ^7 / 



/ ^7 \ 



Fry Hry 0, 



7-1 



C 7 — —Af 1 ^2n — 2A f-y-^h^, 
H=l 

7— 4n 



7-1 



^7 — —Bh~f-2n — A £ e 2az / 7 _4„_ M / M + £ /i 7 _ M /i M , 

/i=0 /^=1 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



(4.11) 
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and 



n 



7 



12A - + a - 1 

n 



2A/ 








V 













7 
n 



-2 -1 



7 



/ 



-12 - 

n 

As expected, the eigenvalues of A 7 are ^ — k r , where k r 
The diagonalizing matrix for A 7 is 



(4.12) 



T, 0, 6, r are the resonances. 



/ " f 1 



(4.13) 



-A/o 1 ^ 

A/ a -A/ (a + 6) 1 - a 

1 3(2a + 5) 

\ -3 12(2a + 5) ) 

with D 1 = P~ 1 A 1 P, and ||P|| ||-P _1 || < M(A, /o), a constant independent of 7, but 
depending upon A and the arbitrary constant /q. Since — ^<A<0or0<A<^, we 
have < f < 5 < 6, so that k r = 6 is again the largest resonance, and all eigenvalues of 
A 7 are positive for ^ > 6. Thus, for ^ > 6, the solution of (4.9) may be expressed as 



/z 
Pe D '< {x - z) p- l F 1 {x)dx, 
-00 



with 



, 7 (z)|| < M(A,/ ) / e(^- 6 )^ 2 )|| J F 7 (x)||^. 



Theorem 4.1. There exists K > suc/i i/tai /or all z < and 7 > 1, 

,i/»i, £ JEL 

As the proof is practically identical to the one of Theorem 3.1, we omit the details. 
By (4.3), 



(4.14) 
(4.15) 
(4.16) 



s = £/ 7 (*)«' 7+an 
7=0 

and by (4.16), 

{my 



\f^(z)w 



7+«n I ^ 



w 



7+on 



and the series 



E 

7=0 



{my 



w 



7+on 



(4.17) 



(4.18) 



(4.19) 



converges absolutely by the ratio test for w < jg, i.e., for r < j^yr- Thus, the series 
(4.17) for x and, similarly, the series (4.5) for y and (4.7) for x and y converge absolutely 
for < r < R, where R is at least j^yi ■ 
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5. Concluding Remarks 

Three types of psi-series solutions of the cubic Henon-Heiles system have been shown to 
be absolutely convergent on intervals of the form < r < R, R > 0. It is clear that 
the series for x, x, y and y converge uniformly on compact subintervals of (0, R), thereby 
justifying the termwise differentiation of the series. 

The resummations that have been performed in order to proceed with the convergence 
proofs amount to the replacement of multiply- indexed series with constant coefficients 
by single series with polynomial coefficients. This procedure makes use of the fact that 
the resonances appearing as exponents within the psi-series have positive real parts. For 
example, to obtain (2.3) and (2.6) from (2.2), it is essential that Re(r) > and Re(f) > 0. 
Otherwise, as occurs in (4.2) when f is pure imaginary, a resummation would give rise to 
series containing f-y(z) which are not polynomials, but infinite series. 

The situation when A = — \ is very interesting. Here, a = —2 (the leading order for 

x), r = 5 and f = 0. But the leading coefficient for x is ±j^2 + j = 0, contradicting the 

fact that x = 0(r~ 2 ). The correct leading-order behaviours have been given in [7] and [9, 
p. 340]. In our notation, 

x ~ (-30) 1/2 r- 2 (ln(r))-5, 

y ~ Zt- 2 + ^t-\Ht))-\ 

A full series expansion of the form 

oo oo 

I = ^aHr"(Mr))-4, 

k=0 1=0 

oo oo 

fc=0 1=0 

fails, due to an incompatible resonance at (k,l) = (0,2). We are not aware of the correct 
series in this case. 
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Appendix A 

This section contains the Lemmas necessary for the convergence proofs. 

Lemma Al. Let 1/0 have constant sign, let q > 1 be an integer, and for 7 > 1, let 

p,{x) = £ c (]>r 

m=0 
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be a sequence of polynomials of degree n 7 = [^]. Given an integer N>q,0<M<l and 
p > 0, there exists K with sign(.K') = sign(X) and p\K\ > 1 such that 

, , Nl (p\K\+ KX)~t/i 

|P 7 (X)|<M ^' M l<7<iV-l. 

This is Lemma 4.4 of [13], to which we refer the reader for a proof. 
Lemma A2. 

7-1 



lim > 

•v — »on ^ — J 



7T. 



The result follows easily by the proof of the integral test for convergence of series. For the 
details, we refer to Hemmi [19, p. 64]. 

Lemma A3. Let X, Y > and for 7 > 1, let P 1 (X,Y) be a sequence of polynomials of 
degree 7 in X and Y, i.e., 

P 7 (X,Y) = j2Y.c$X»Y^. 

/3=0/i=0 

Given an integer N>l,0<M< \ and p > 0, there exists K > such that 

IW)I < M <<™™: for 1 < 7 <^_ 1 . 

Proof. Let 

1/7 



.ftT = max 

1 < 7 < at- 1 

< (3 < 7 
< /i < /3 



7 \(/? 



Mp-t-f 3 



then for 1 < 7 < iV - 1, < /? < 7, < n < /3, we have 
(7), , M ( l\( P\ r^~,-0 



which implies that 

7 

\P y (X,Y)\<^2^T\c^\X»Y^ 

(3=0 n=0 

7 P 



< yy( i)( P) r^-W-" 



M 



/3=0 M=0 

(pit - + kx + irr) 7 



by the binomial theorem. 
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Appendix B 

First, we show that the solution (2.1), namely, 



oo oo oo 

-EEE«'^" rl+rm 

k=0 1=0 rn=0 



oo oo oo 



(Bl) 



» = EEE w*" 2+ '''"''"' 

k=0 1=0 m=0 



of the system (1.4), corresponding to the case(i)-leading orders and non-integral reso- 
nances, is self-consistent and contains four arbitrary constants. 
Substitution of (Bl) into (1.4) gives 



3/1 3 
a 00 o = ±^y2 + -, 6 000 = --, (B2) 

and for (k,l,m) / (0,0,0), the coefficient recursion relations 
(k + rl + fm){k + rl + fm - 5)a k i m + 2Aaooo&fcZm 



— Aa>k—2,l,m 2A ^ ^ 0,pqsbk—p,l—q,m—sj 

(0,0,0)<(p,q,s)<(k,l,m) 



2Xa 000 aki m + 



6 

(k + rl + rm — 2)(k + rl + rm — 3) + — 

A 



(B3) 



— Bbk—2,l,m ^ ] [ ^dpqsQ'k— p,l— q,m— s ~l~ ^pqs^k— p,l— q,m— s] j 

(0,0,0)^(p,q,s)x(fe,;,m) 

where we have employed the notation (0, 0, 0) -< (p, g, s) -< (k, I, m) to mean that < p < k, 
0<q<l,0<s<m, with (p, q, s) / (0, 0, 0) and (p, q, s) / (k, I, m). Expressed in matrix 
form, the left-hand sides become 



(k + rl + rm)(k + rl + rm — 5) 2Aaooo 

2Aa 00 o (k + rl + fm — 2){k + rl + fm - 3) + 6 



O'klm 
him 



,(B4) 

A 

and the coefficient matrix is singular precisely when k + rl + fm has one of the values 
— 1, 6, r, f, i.e., when 

(fc, I, m) = (-1, 0, 0), (6, 0, 0), (0, 1, 0), (0, 0, 1). (B5) 

These are the resonances of (Bl), and compatibility must be checked at these values, 
confirming both the self-consistency of (Bl) and the presence of four arbitrary constants. 

The resonance at (k,l,m) = (—1,0,0) is trivially compatible, and corresponds to the 
arbitrariness of to- At (k, l,m) = (0, 1, 0), we find that 

a io is arbitrary, 6 io = — «oio, (B6) 

^Aaooo 

and at (k, I, m) = (0, 0, 1), that 

y I y 5 J 

a oi is arbitrary, 6 oi = — «ooi- (B7) 

^Aaooo 
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At (k,l,m) = (6,0,0), the compatibility condition is 

Aaooo040o^. + 2A 2 0000(^200^400 + 04006200) — 36400-B — 6(Aa2ooa40o — 62006400) = 0.(B8) 
The coefficients which affect (B8) are 



1 3 

oooo — +^V/ 2 + 0000 — — ^, 

OlOO = ^300 = ^500 = 0, 6100 = 6300 = 6500 = 0, 

aooo (x + B) B - A (2 + ±) 

G20 ° = 12(1+1) ' b20 ° = 4A(1 + I) ' (B9) 
Q200 (l + j) (A + 2A6 2 qq) - Aapoo (-66200 + Aa 2 ^ - b|oo) 

10(4+1) ' 

Aaoooa20o( j 4 + 2A6200) — 2 (-B6200 + Aa^oo ~~ ^200) 



«400 



(BIO) 



10 (4 + 3 ) 

from which (B8) is confirmed. The system governing a6oo an d 6600 is 

(6 2Aaooo \ / a 600 \ / — ^ a 400 — 2A(a2oo^40o + 04006200) 
2Aa 000 12 + f / V 6 eoo J \ - B b mo - 2(\a 200 a mo - 62006400) 
and reduces to 

6a600 + 2Aaooo660o = —^400 — 2A(a2oo640o + 04006200)1 (Bll) 

showing that agoo is arbitrary, and 6600 is defined in terms of agoo- 

Next, we demonstrate that (2.1) is equivalent to (2.2). The series which defines x in 
(2.1) may be broken up into three parts: m = I, m < I, and m > I. Thus, 

00 00 00 00 i— 1 

* = EE «M;r fc - 2+5 ' + E E E «^- 2+5m+( '- m)r 

k=0 1=0 k=0 1=1 m=0 

(B12) 

00 00 m— 1 
fc=0m=l 1=0 

where use has been made of the fact that r + f = 5. The first sum in (B12) is expressible 
as 

00 00 00 00 

E E ^ur k ~ 2+51 = E E a ^~ 2 = E «^~ 2 , (B13) 

k=0 1=0 i=0 k+5l=i i=0 



where 



aw = ^2 a kii- (Bl 4 ) 

k+hl=i 
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The second sum in (B12) may be manipulated as follows: 

oo oo I— 1 

EEE ^l m r k ' 2+5m ^ l ~ m > 



k=0 1=1 m=0 

oo oo / 



EEE^mV^^ (luting j = l-m) 

k=0 1=1 j=l 



oo oo oo 



y~] ^ cik : i : i-jT k 2+5(7 fi + i r (reversing the order of summation) 

k=0 j=l l=j 



oo oo oo 



= E E E ak,l +j ,lT k - 2+5l+jr (replacing I by I + j) 
fc=o j=i z=o 

oo oo 

= EE E 

i=0 j=l k+5l=i 



= E E a ^ r 



2+jr 
i=0 j=l 



where 



(B15) 



a u = E a M+j,«- ( B16 ) 

fc+5/=i 

Similarly, the third term in (B12) becomes 

oo oo 

EE^' f - ( B17 ) 

i=0 j=l 

where 

«ij = E a k,m,m+j- (B18) 
fc+5m=i 

Combining (B13), (B15) and (B17), we obtain 

oo oo oo oo 

* = E E «^~ 2+jr + E E «^'- 2+J '- (Bi9) 

i=o i=o i=o j=i 

Similarly, 

oo oo oo oo 

» = E E h ^- 2+jr + E E ^ l_2+Jf • ( B2 °) 

i=0 j=0 i=0 j=l 

Defining a i0 = b i0 = for all i > 0, (B19), (B20) coincide with (2.2). 
Finally, the solution given in [3] is 

oo oo oo oo 

x = E ^- 2+3{r - 2) + E E ^i^" 2+j(r " _2) (B2i) 

j=0 i=0 j=l i=0 
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(their equation (3.5a), in our notation), with a similar expression for y. To obtain the 
form (B21) from (B19), write 

oo oo oo oo 

i=0 j=0 i=0 j=l 

oo oo oo oo 



E E a^ +2j - 2+j(r - 2) + E E out* 2 *- 2 **?- 2 ) 

i=0 j=0 i=0 j=l 

oo oo oo oo 

E E ^r 1 - 2 ^ + E E Si-^r*-^- 2 ) 

j=0 i=2j j=l i=2j 

oo oo oo oo 

E E Atf^ 7 **-** + E E ^■ir i - 2+j ' (f - 2) , 

jr=0 j=0 jr' = l 1=0 



(B22) 



where 



A 3l = \ «-™> if J ^ 1 , ^, = ( ''' 2 "' * * *V \, (B23) 
' 0, otherwise J ' J \ 0, otherwise ' 



and similarly for y. 

Appendix C 

First, we confirm that the solution (4.2), namely, 



oo oo oo 



^ = EEE a ^ fe+s+(2+s)m+rt ' 



fc=0 £=0 m=0 

oo oo oo 



(CI) 



!/ = EEE^" 2+(2+s)m+f ' 



fc=0 1=0 m=0 

of the system (1.4), corresponding to the case(ii)-singularities, is self-consistent and con- 
tains four arbitrary constants. 

Substitution of (CI) into (1.4) gives 

a oo is arbitrary, 6 oo = 6, (C2) 

and for (k,l,m) / (0,0,0), the coefficient recursion relations 

[k + (2 + a)m + fl][k + (2 + a)m + r(l - l)]a klm + 2Aa oo^m 

= — ^4 a fc— 2,l,m ~ 2A ^ ^ a k—p,l—q,m—rbpqri 

(0,0,0)-<(p,q,r)-<(k,l,m) 



[k + {2 + a)m + fl + lp + (2 + a)m + fZ - 6]6 fcZm 

= —Bbk-2,l,m — A ak-p,l-q, m -r-20'pqr 
p=0 q=0 r=0 

~\~ ^ ] bk—pj—q^m—rbpqr-, 

(0,0,0)-<(p,q,r)-<(k,l,m) 



(C3) 
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where we have defined the symbol as in Appendix B. Expressed in matrix form, the 
left-hand sides become 



[k + (2 + a)m + fl] [k + (2 + a)m + f(l - 1)] 2A(i000 

[k + (2 + a)m + fl + 1] [fe + (2 + a)m + FZ - 6] 



(C4) 



and the coefficient matrix is singular precisely when k + (2 + a)m + W has one of the 
values — 1,6, f, i.e., when 

I, m) = (-1, 0, 0), (6, 0, 0), (0, 1, 0). (C5) 

These, together with (k,l,m) = (0,0,0), are the resonances of (CI). Compatibility at 
(0, 0, 0) and the arbitrariness of aooo has already been confirmed. The compatibility of 
(—1,0,0), reflecting the arbitrariness of to, is trivial, as usual. It remains to check the 
remaining two. 

At (k,l,m) = (0, 1,0), we find that 

^oio = 0, aoio is arbitrary. (C6) 

The compatibility condition at (k, I, m) = (6, 0, 0) is 

5 

- -B&400 + ^2 h-pflfibpoo = 0. (C7) 
P =i 

The coefficients which affect (C7) are 

B B 2 

6100 = 0, ^300 = 0, ^500 = 0, ^200 = -^-> ^ 400 = 4o"' 

from which (C7) is confirmed. The second equation in (C3) then shows that ^600 is 
arbitrary, and the first gives a6oo hi terms of beoo. 

Second, we note that a solution different from (CI) has been given in [3], but we find 
that the coefficient recursion relations (their equation (3.9)) do not follow from the solution 
(their equation (3.6)) without neglecting certain terms, precisely those terms which lead 
to anomalies when a resummation is performed in an attempt to prove convergence. 

Finally, we remark upon the "derivation" of the form of (CI). The exponent fl is 
required in order to produce an arbitrary coefficient corresponding to the resonance at f. 
This is usual. The exponent am is required in order that all terms in the system (1.4) can 
be balanced. But since Re(a) > —2, am must be replaced by (2 + a)m in order to avoid 
exponents of r with arbitrarily large negative real parts, which would render the series 
invalid. 
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